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ABSTRACT 

The observed redshift distribution of Ly-a lines and metal systems is examined in 
order to discriminate and to trace the evolution of structure elements observed in 
the galaxy distribution, at small redshifts, and to test the theoretical description of 
structure evolution. We show that the expected evolution of filamentary component of 
structure describes quite well the redshift distribution of metal systems and stronger 
Ly-a lines with log(Njji) >14, at z < 3. The redshift distribution of weaker Ly-a lines 
can be attributed to the population of poorer structure elements (Zel'dovich pancakes), 
which were formed at high redshifts from the invisible DM and non luminous baryonic 
matter, and at lower redshifts they mainly merged and dispersed. 



Key words: cosmology: large-scale structure of the Universe 
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1 INTRODUCTION 

During the last years an essential progress has been achieved 
both in the observations and interpretations of absorption 
spectra of quasars. Now a representative set of absorption 
lines with measured redshifts, column densities, and Doppler 
parameters is available and the main characteristics of ab- 
sorbers and their evolution are, in general, established (see, 
e.g., Carswell 1995; Cristiani et al. 1995, 1996; Hu et al. 1995; 
Fernandez-Soto et al. 1996; Ulmer 1996; Cooke et al 1997; 
Dinshaw et al. 1995; Bahcall et al. 1993, 1996; Kim et al. 
1997). The available information suggests various descrip- 
tions of absorbers ranging from pressure-confined gaseous 
clouds to associations with halos of galaxies (see discussions 
in Rees 1995; Charltone 1995; Fernandes-Soto et al. 1996; 
Miralda-Escude et al. 1996; Petitjean 1997). If at low red- 
shifts a significant number of stronger Ly-a lines and metal 
systems are certainly associated with galaxies (Bergeron et 
al. 1992; Lanzetta et al. 1995; Cowie et al. 1995; Tytler et al. 
1995; Le Brune et al. 1996) than the population of weaker 
absorbers, dominating at higher redshifts, has practically 
disappeared by the redshift z « 2 and only some weak Ly-a 
lines are observed far from galaxies even at small redshifts 
(Morris et al. 1993; Shull 1997). 

Now the origin and evolution of absorbers are linked to 
the process of formation and evolution of high density struc- 
ture elements (Petitjean et al. 1995; Miralda-Escude et al. 
1996; Hernquist et al. 1996; Zhang et al. 1997, 1998; The- 
uns et al. 1998, 1999; Bryan et al. 1999; Dave et al. 1999; 



Weinberg et al. 1998; Machacek et al. 2000). The numerical 
simulations of dynamical and thermal evolution of gaseous 
component use the same spectrum of initial perturbations 
that is responsible for the formation of DM and observed 
galaxy spatial distributions, and they successfully reproduce 
the main observed properties of absorbers. This allows us to 
consider the evolution of absorbers in a context of a more 
general process of nonlinear evolution of small initial pertur- 
bations of DM component and formation of observed galaxy 
distribution and to compare theoretical conclusions with ob- 
servations. 

This approach links the properties of absorbers, at high 
redshifts, with the spatial galaxy distribution at small red- 
shifts and allows to trace some observational characteristics 
of nonlinear evolution of matter distribution up to redshifts 
z ~3. 

The analysis of large modern redshift surveys such as 
the Durham/UKST Galaxy Redshift Survey (Ratcliffe et al., 
1996) and the Las Campanas Redshift Survey (Shectman et 
al., 1996) demonstrates that the galaxy distribution in the 
universe can be roughly described as a joint network com- 
posed of galaxy filaments and walls. The analysis of both 
surveys shows that about 50% of galaxies are concentrated in 
wall-like elements surrounding huge underdense regions with 
a typical size of 50 - 100ft" 1 Mpc (h = H o /100 km/s/Mpc 
is the dimensionless Hubble parameter) (Doroshkevich et al. 
1996, hereafter LCRS1; Doroshkevich et al. 2000). In many 
respects such elements are similar to the Great Wall (de Lap- 
parent, Geller & Huchra 1988; Ramella, Geller, & Huchra 



© 0000 RAS 



2 Demianski et al. 



1992). Galaxy filaments appear between walls and intertwine 
them into a joint structure. The hierarchy of poorer and 
sparser filaments - including single galaxies - can be con- 
sidered as a bridge between the richer and poorer structure 
elements formed by the non luminous baryonic and dark 
matter (DM). Such poorer elements contain some fraction 
of hot gas and could be seen as Ly-a absorbers away from 
galaxies. Observations of weak Ly-a absorbers in voids, men- 
tioned above, can be attributed to those poorer DM struc- 
ture elements. 

The evolution of structure elements formed by dark 
matter was investigated both theoretically, using the 
Zel'dovich theory (e.g., Zel'dovich 1970; Zel'dovich & 
Novikov 1983; Shandarin & Zel'dovich 1989), and numer- 
ically (e.g., Shandarin et al. 1995; Cole et al. 1997, 1998; 
Jenkins et al. 1998; Doroshkevich et al. 1999, hereafter 
DMRT; Demianski et al. 2000, hereafter DDMT). The sta- 
tistical description of DM structure evolution for the CDM- 
like power spectra was given in Demianski & Doroshkevich 
(1999 hereafter DD99) and in DDMT. The main simulated 
characteristics of walls are found to be consistent with the- 
oretical expectations. 

Some progress in the theoretical interpretation of prop- 
erties of absorbers has been recently achieved by Hui et al. 
(1997) and Nath (1997), who used the Zel'dovich approxi- 
mation for a semi analytical description of some of the ab- 
sorbers characteristics. Another considered theoretical ap- 
proaches were based on the 'minihalo' model (Meiksin 1994), 
'Cosmic Webb' theory (Bond & Wadsley 1997) and Press- 
Schechter method (Valageas et al. 1999). 

In this paper we examine the observed redshift distribu- 
tion of absorbers using the available information about the 
structure parameters at small redshifts (LCRS1; Doroshke- 
vich et al. 2000; DMRT) and statistical description of DM 
structure evolution (DD99, DDMT) in models with CDM- 
like power spectra. We assume that the Lyman-a clouds 
trace the potential wells formed by DM pancakes, filaments 
and walls. This assumption was widely discussed earlier 
(Rees 1986; Ikeuchi & Ostriker 1986; Bond, Szalay & Silk 
1988; Miralda-Escude & Rees 1993; Meiksin 1994; Hui et 
al. 1997; Nath 1997; Bond & Wadsley 1997; Valageas et al. 
1999), and it is qualitatively consistent with results of nu- 
merical simulations (Petitjean et al. 1995; Hernquist et al. 
1996; Petitjean 1997; Theuns et al. 1997; Dave et al. 1999; 
Weinberg et al. 1998). 

The potential of this approach is limited due to the 
small number of observed absorbers and because the char- 
acteristics of structure derived from observations are not 
sufficiently accurately known. The theoretical description of 
DM structure evolution is also not complete. In particular, it 
cannot yet describe adequately the important process of dis- 
ruption of structure elements and their transformation into 
a system of high density clumps. This process depends on 
many factors (Doroshkevich 1980; Vishniac 1983; Valinia et 
al. 1997) and strongly influences the observed properties of 
absorbers both at large and small redshifts (Miralda-Escude 
et al. 1996; DMRT; DDMT). Moreover, the evolution of DM 
structure elements and observed absorbers is not identical. 

Nonetheless, even so limited information allows us to 
obtain a reasonable description and interpretation of some 
of the observed characteristics of absorbers. In this paper 
we consider the statistical description of redshift distribu- 



tion of absorbers and metal systems which is more detailed, 
in some respects, then that usually used. The comparison 
of measured redshift distributions with approximate theo- 
retical expectations (DD99) allows us to discriminate sta- 
tistically the contributions of high density filamentary and 
lower density sheet-like structure elements formed by DM 
and gaseous components of matter, and to trace their evo- 
lution. These two components of structure elements can be 
roughly identified with two subpopulations of absorbers dis- 
cussed, for example, in Bahcall et al. (1996) and Weymann 
et al. (1998). 

This division of structure elements into filaments and 
more sheet-like walls and pancakes is inevitably statistical 
only. The more detailed analysis of DM and galaxy distribu- 
tion, at small redshifts (DMRT, DDMT), demonstrates lim- 
itations of such description and shows that usually the mor- 
phology of structure elements can be more accurately char- 
acterized in terms of degree of filamentarity and/or sheetness 
only. The investigation of morphology and spatial distribu- 
tion of DM and gaseous structure elements in simulations at 
high redshifts with the new powerful techniques (such as the 
Minimal Spanning Tree or Minkowski Functional), as well as 
direct comparison of morphological and other characteristics 
of the structure elements with their redshift distribution is 
required to clarify this problem and to obtain more detailed 
comparison of observed and simulated matter distributions. 

This paper is organized as follows. The theoretical 
model of structure evolution and the technique necessary 
for the analysis are briefly described in Sec. 2. Sec. 3 con- 
tains information about the used observational databases. 
The results of statistical analysis are presented in Sec. 4. 
Discussion and conclusion can be found in Sec. 5. 



2 REDSHIFT DISTRIBUTION OF 
ABSORBERS 

2.1 Pancakes, filaments and walls associated with 
absorbers 

The redshift distribution of absorbers can be characterized 
by the mean comoving free path between absorbers < / > 
(Sargent et al. 1980) as 

dN(z) =< I y- 1 c(l + z)dt =< I >~ x ^4 (2.1) 

H(zj 

H(z) = H ^n m (l + zf + (1 - n m - n A )(l + z) 2 + Ov 

where dN(z) is the number of absorbers between z and 
z + dz, < I > _1 is the mean comoving linear number den- 
sity of absorbers, c is the speed of light, Q m and Ov arc 
correspondingly the dimensionless matter density and cos- 
mological constant. The redshift dependence of the mean 
comoving linear number density of absorbers can be written 
as: 

n a bs = -£r- < I > _1 = n w (z) + rif(z) + n pan (z), (2.2) 
no 

dN(z) = n *bs-£jj^, 

where the three terms in (2.2) correspond to three main 
types of structure elements observed at small redshifts: the 
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rich walls (n w ), filaments (n/), formed by galaxies, and poor 
pancakes (n pan ), situated far from observed galaxies. 

Such division of structure elements into three different 
subpopulations is made for the sake of simplicity only. In 
particular, there is an anisotropic sheet-like low density halo 
around all filaments, which in turn are later strongly dis- 
rupted into a system of high density clumps. The higher 
density ridges cross often the low mass pancakes. Walls arc 
also often formed by a system of high density clumps and 
filaments connected by lower density bridges. This means 
that in nature, there is really only one population of struc- 
ture elements with a broad variety of continually distributed 
properties, and these elements can be considered as the more 
typical and conspicuous ones. Nonetheless, the considered 
classification makes this analysis simpler and more trans- 
parent, as the rate of structure evolution depends upon the 
degree of filamentarity or sheetness. 

2.2 Evolution of DM structure elements 

In the DM dominated universe the basic characteristics of 
matter distribution are determined by the process of forma- 
tion and evolution of DM structure elements. The observed 
characteristics of absorbers and parameters of DM structure 
elements are not identical but they are closely connected and 
relations between these parameters can emerge from suitable 
investigations. 

For the CDM-like transfer function (Bardeen et al. 
1986) and the Harrison-Zel'dovich initial power spectrum 
the approximate statistical description of the formation and 
evolution of DM structure elements based on the statistics 
of initial perturbations was obtained in DD99. The random 
formation and merging of Zel'dovich pancakes with 'mass' 
(column density of a pancake) m and their evolution, due to 
the transversal expansion or compression, are described by 
approximate expressions similar to the Press-Schechter rela- 
tions. Similar approximate relations can be also written for 
the filamentary component of the structure. Here we will ap- 
ply some of these results to describe the observed evolution 
of absorbers. 

As is usual in the Zel'dovich theory the redshift depen- 
dence of all characteristics of structure is expressed through 
a function B(fl m ,z), which describes the growth of per- 
turbations in the linear theory. For the flat universe with 
fi m + ft a = 1, fim > 0.1, the function B(z) can be ap- 
proximated (DD99) with a precision better then 10% by the 
expression 



B(z) 



_ 3 1 -O m + 2.2^(1 + 2 ) 3 



1 + 1.2fi m 

and for an open universe with 0.1 < f2 m < 1, Oa = as 



(2.3a) 



D/ \ — 1 1 1 2.5Qm 



(2.36) 



(Zel'dovich & Novikov 1983). For Q m = 1, fi A = both 
expressions give B^ 1 (z) = 1 + z. 



2.2.1 Formation and evolution of DM structure 

The evolution of DM structure can be outlined as a progres- 
sive matter concentration within more and more massive 
structure elements which can be roughly discriminated into 



pancakes, filaments and walls. The structure elements are 
primarily formed as Zel'dovich pancakes (Shandarin et al. 
1995), but later are successively transformed into filaments 
and/or (elliptical) clouds. At high redshifts 2 > 5 - 7 a sig- 
nificant fraction of matter (~80%) can be compressed into 
low mass pancakes with M ci ~ 10 7 - 10 9 M Q . Later some of 
such clouds can evolve into dwarf galaxies, some could dis- 
perse, but the main fraction will be accumulated by more 
massive structure elements. Properties of the surviving and 
relaxed clouds can be similar to those discussed in the 'mini- 
halo' model (Rees 1986, 1995; Miralda-Escude & Rees 1993; 
Meiksin 1994) 

The filamentary component of structure also forms at 
high redshifts due to the compression of some part of pan- 
cakes, and, at z ~ 1 - 3, it is, probably, the most conspicuous 
population of high overdensity structure elements accumu- 
lating 40 - 50% of DM particles (see., e.g., Governato et al. 
1998, Jenkins et al. 1998). Further evolution of structure can 
be described as a successive merging of lower mass pancakes 
and filaments into more and more massive wall-like structure 
elements. The largest observed wall-like elements, similar to 
the Great Wall, are formed at redshifts z <1 due to the infall 
and intermixture of smaller pancakes and filaments. 

In observations, at 2 < 1, galaxy filaments are usually 
situated within extended under dense regions surrounded 
by higher density walls, and they form an irregular bro- 
ken network. This network can be described as a trunk, 
that is the longest filament of the complex, and a set of 
shorter branches. This means that the spatial distribution 
of filaments can be conveniently characterized by the mean 
2D surface number density, o~f(z), that is the mean num- 
ber of filaments intersecting a unit area of arbitrary ori- 
entation. For observed galaxy filaments it is estimated as 
(7/(0) « (1 - 0.7) ■ 10- 2 fe 2 Mpc~ 2 (LCRS1). The spatial dis- 
tribution of observed walls can be characterized by their 
mean separation (D w ) and by their linear density, a w , that 
is the mean number of walls intersecting a unit length of 
a random straight line. For the observed galaxy walls it is 
estimated as (D w ) «40 - 60fe _1 Mpc that corresponds to 
a w Ri2.5 - 1.6-10 _2 /iMpc _1 . The low mass pancakes can be 
considered as isolated objects and characterized by the mean 
3D number density v van {z). 

Simulations show that rich pancakes and filaments are 
relaxed (at least along shorter axes) and are surrounded 
by an extended DM (and gaseous) halo. Central parts of 
dense elements are gravitationally confined and are unstable 
with respect to the small scale disruption (DDMT) . Some of 
the high density clumps are usually embedded in filaments 
and walls what essentially accelerates the relaxation of com- 
pressed matter. In observations, effects of such disruptions 
are seen as groups and clusters of galaxies. 

The basic properties of DM structure elements vary in a 
broad range, but it can be expected that galaxies are formed 
mainly within richer high density filaments and pancakes. 
The essential fraction of poor pancakes and, perhaps, fila- 
ments, can be invisible, as it does not contain sufficiently 
bright galaxies. Such structure elements can be observed as 
weak Ly-a absorbers situated far from any galaxy. Some 
number of such absorbers situated at a distance of up to 
5h~ Mpc from the closest galaxy was found by Morris et al. 
(1993) and Shull (1997). 

This description agrees well with results of numerical 
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simulations at small redshifts (DDMT). The validity of this 
description, when it is applied to structure at high redshifts, 
is not yet reliably verified with available simulations due to 
a small density contrast of poor structure elements. The first 
tests show, however, the existence of three kinds of structure 
elements, namely, high density filaments and clumps, and 
low density pancakes in DM spatial distribution at z =3. 
The theoretical description (DD99) confirms also the self- 
similar character of structure evolution (at least when the 
Zel'dovich approximation can be applied). This picture is 
also qualitatively consistent with other simulations of struc- 
ture evolution (Miralda-Escude et al. 1996; Governato et al. 
1998; Jenkins et al. 1998; Meiksin 1998; DMRT; Dave et al. 
1999; Weinberg et al. 1998). 

2.2.2 Main factors of structure evolution 

The four main factors characterizing the structure evolution 
are 

(i) successive formation of structure elements with larger 
and larger masses, 

(ii) merging of earlier formed, more numerous, low mass 
filaments and pancakes 

(iii) expansion or compression of matter along filaments 
and pancakes 

(iv) small scale clustering of compressed matter and dis- 
ruption of structure elements. 

The action of these factors was discussed in DD99 (Sees. 4.1, 
4.2 & 5.2) and it was shown that for the CDM-like initial 
power spectrum and Gaussian initial perturbations it can be 
approximately described by simple analytical expressions. 
Some of these relations were tested and confirmed through 
the direct comparison with DM simulations (DDMT). Here 
we will briefly summarize these results. 

The distribution function of Zel'dovich pancakes with 
mass (surface density) m p , N p (m p ), can be written as fol- 
lows: 



N p (m p ) oc 



\J m p {m p ) 



exp 



(m p ) 



erf 



(m P ) 



, (2-4) 



(m p ) oc B , 

where the successive merging of low mass pancakes in the 
course of formation of richer pancakes is described by the 
er /-function. The process of merging decreases the number 
density of numerous low mass structure elements - both fil- 
aments and pancakes - with m p <C (m p ), but its influence 
becomes negligible for rare, massive, wall-like structure el- 
ements, with m p > (m p ), and rare rich filaments. For the 
fraction of pancakes with m m i n < rn p < m max from (2.4) 
we have: 



fv « er f 



erf 



(m P ) 



(2.5) 



For rich walls formed presumably at z < 1 — 1.5 with 
mmax — ► oo and m min = m w ~ (m p ) this relation is simpli- 
fied and approximately we have 



f w oc 1 - erf 





/ m w 


oc erfc 


1 m w 


A 


1 (m p )_ 


y (m p )_ 



(2.6) 



what constrains the population of massive pancakes. For the 
subpopulation of poor pancakes similar relation 



f w oc 1 - erf 



(m P ) 



oc exp 



(m p ) 



is valid only at higher redshifts z > 3.5 - 4, when m max S> 
m m in ~ (m p ), whereas during the most interesting period, 
z < 3.5, when m min <C (m p ), the second term in (2.5) can 
be omitted. 

The discrimination of rich walls and poor pancakes 
seems to be artificial, in some respects, but it allows us to 
provide a closer connection with observed galaxy and simu- 
lated DM distributions at smaller redshifts. Moreover, rich 
walls are certainly strongly nonhomogeneous and disrupted 
into high density clumps and filaments, what essentially de- 
creases their covering factor. Due to these peculiarities such 
walls can form a special class of objects. 

The distribution function of filaments with mass (linear 
density) m/ is found to be proportional to 



N f oc 



1 



(m f ) 



m f \ 
4{m f ))> 



(mf) oc B 4 



(2.7) 



where Ko(x) is the modified Bessel function (DD99) and the 
fraction of rich filaments, with m./ > (m/), is 



// oc exp 



m f 



(m/) J 



(2.8) 



As before, the infall of filaments into rich walls decreases 
the fr action of filaments associated with the network oc 
erf[^m p /{m p )] oc B' 1 . 

Furthermore, the properties of structure are essentially 
changed due to matter expansion and/or compression in 
transversal directions. Compression results in transforma- 
tion of pancakes into filaments and high density clouds, what 
decreases their effective surface and, also, the probability 
to see such a structure element as an absorber. Strong ex- 
pansion increases the surface of pancakes and the length 
of filaments and correspondingly decreases their density. In 
the limiting case, low mass pancakes and filaments can be 
transformed into a system of high density clumps or they 
can even completely disperse. The expansion decreases also 
the column density of HI below the observational threshold, 



N 



(thr) 



what makes such elements invisible as 



absorbers. Therefore, the population of structure elements 
identified with absorbers is restricted by the condition of 
their relatively slow expansion or compression in transver- 
sal directions. Fract ion of suc h low mass pancakes decreases 
with time oc erf 2 [y/ m p / (m p )], and fraction of such low mass 
filaments decreases oc erf[yj m p /{m p )] oc B^ 1 . The influ- 
ence of these factors is less important for rare rich walls and 
filaments. 

Due to gravitational instability of compressed matter 
both pancakes and filaments are usually disrupted into a 
system of high density clouds linked by low density bridges. 
Such disruption occurs more rapidly in the central high den- 
sity parts of structure elements, but, probably, its impact is 
not so important for the extended lower density halo. When 
these processes are taken into account the final relations 
become very cumbersome and therefore such more detailed 
treatment is not really justified for when only limited set of 
observed absorbers is available. 
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When these factors are taken into account we can only 
approximately describe the expected evolution of walls, fila- 
ments, and low mass pancakes associated with the absorbers. 
The comoving linear number density of massive walls, a w (m) 
with m > m w , the comoving surface number density of fil- 
aments, <r/(m), and the comoving number density of pan- 
cakes, v P an(m), can be approximated by simple expressions: 



a w oc exp 



Of cc B exp 



b w {m) 
B 2 



bf(m) 
B' 2 



(2.9) 



(2.10) 



(2.11) 



where the parameters b w ,bf k, b^m) and the exponential 
terms in (2.9) and (2.10) restrict formation of structure el- 
ements with masses > m, the term B~ 2 in (2.10) describes 
successive infall of filaments into richer walls and the de- 
crease, due to expansion or compression, of number of fila- 
ments observed as absorbers with Nhi > N^ r \ The de- 
crease of the number density of pancakes is described by the 
factor bi in (2.10). 

For restricted redshift interval z < 3.5 the three param- 
eter expression for v van can be approximated by a simpler 
one parameter expression 



oc B exp — 



*(m) 



B 2 



(2.12) 



which correctly describes the asymptotical behavior of v vari 
at small and large redshifts. In the intermediate region it 
provides a reasonable precision ~ 10 - 15%, for z < 3.5. The 
more complicated relation (2.11) provides better description 
of absorbers evolution especially at higher redshifts z > 3 - 
3.5 and can be used for more refined fits. 

High density clouds and galaxies are usually embedded 
in filaments and pancakes and cannot be distinguished as 
a special class of absorbers. These problems were discussed 
with more details in our other publications DD99, DMRT 
and DDMT. 



2.3 Model of absorbers evolution 

The relations (2.3), (2.9) - (2.12) cannot give the full de- 
scription of observed redshift distribution of absorbers, but 
they suggest possible redshift evolution of various types of 
absorbers and introduce the function B(z, f2 m ) as an impor- 
tant characteristic of such evolution. It is important that the 
redshift evolution of free-path between walls, filaments, and 
low mass pancakes is expected to be different, and there- 
fore, it might be possible to single out statistically these 
three subpopulations of absorbers, and to establish correla- 
tions between the evolutionary rates and other observational 
characteristics of absorbers. 

The DM structure is quite complicated in itself as it 
is composed of several types of structure elements with dif- 
ferent evolutionary histories. Furthermore, properties and 
evolution of the observed gaseous and DM structure ele- 
ments are not identical, because of the influence of addi- 
tional factors on the evolution of gaseous component of DM 



confined structure elements and, therefore, even more com- 
plicated evolution of absorbers can be expected. The action 
of these factors usually decelerates the evolution of gaseous 
component but, due to very small density of this component, 
its impact on the evolution of DM component is negligible. 
Thus, the formation of weak absorbers within "minivoids" 
does not accompanied by the formation any DM structure 
elements (Bi & Davidsen 1997; Zhang et al. 1998; Dave et 
al. 1999). 

The observational restrictions such as the condition 



> N 



(thr) 



constrain also the population of absorbers 



and the observed evolution of absorbers. Thus, for example, 
the exponential growth of fraction of massive high tempera- 
ture pancakes (in the limiting case, walls similar to the Great 
Wall) is not found in the observed distribution of absorbers 
(Sec. 4). Moreover, the mean measured b - parameter of ab- 
sorbers, related to the gas temperature, remains constant for 
the observed redshifts 2.5 < z <3.3. This means that, possi- 
bly, some rich DM structure elements are not yet identified 
in the absorption spectra of quasars. 

The redshift dependence of n a b s is driven both by the 
evolution of spatial characteristics of DM structure elements, 
discussed above, and by evolution of their covering factor, 
which characterizes the probability of formation of absorp- 
tion line at the intersection of line of sight and a structure 
clement. The covering factor depends on local conditions 
such as, for example, variations of the UV background and 
activity of nearest galaxies. Now it can be also described 
phenomenologically. 

The important role of UV background is well estab- 
lished after the discussion in Sargent et al. (1980), but 
the UV intensity and its possible variation with redshift 
is known only with large uncertainty (see, e.g., Haardt & 
Madau 1996; Cook, Espey & Carswell 1997). The influence 
of systematic variations of UV radiation is more important 
at small (z < 2) and higher (z > 3 - 3.5) redshifts where 
the statistic of absorbers is strongly limited. At intermediate 
redshifts 2.5 < z < 3.5 where the main fraction of observed 
absorbers is situated the systematic variations of UV radi- 
ation are not so strong, but this radiation is probably re- 
sponsible for the significant irregular variations of absorbers 
density. 

The redshift distribution of observed absorbers can be 
distorted by the formation of artificial caustics in the red- 
shift space (McGill 1990; Levshakov & Kegel 1996, 1997). As 
was shown in DDMT the impact of this effect is certainly 
small at small redshifts. At intermediate redshifts and for 
the usually used CDM-like power spectra the formation of 
artificial caustics is partly suppressed due to strong matter 
concentration within low mass structure elements discussed 
in DD99 and Zhang et al. (1998). This factor transforms the 
continuous matter infall into pancakes to discontinuous one, 
increases the density gradient near pancake boundaries and 
partly prevents the formation of artificial caustics. Moreover, 
such an artificial caustic is actually the preliminary stage of 
formation of a real caustic, and, so, it is rapidly transformed 
into a real one. This factor can moderately change the dis- 
cussed redshift dependence of absorbers. 

The contribution of artificial caustics and absorbers 
within "minivoids" as well as of short-lived rapidly expand- 
ing or contracting pancakes is more important at higher red- 
shifts, when weak absorbers dominate. These factors gener- 
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ate an essential noise what is usually typical for the period 
preceding the epoch of regular evolution. At such redshifts, 
the available statistics of observed absorbers is very poor, 
and our analysis becomes unreliable. Perhaps, more detailed 
and careful investigations based upon the simulations will al- 
low to discriminate observationally between this noise and 
long-lived structure elements (see, e.g., discussion in Zhang 
et al. 1998) that will essentially improve and simplify the 
comparison with theoretical expectations. 

The more detailed description of structure evolution sig- 
nificantly increases the number of fit parameters for limited 
available statistics of observed absorbers, and more detailed 
description accompanied by a further growth of number of 
fit parameters does not seem to be justified. Hence, as the 
first step, we will fit the observed redshift distribution of 
absorbers to relations similar to theoretically expected ex- 
pressions (2.9)-(2.12), but instead of the functions b(m), 
arbitrary fitting parameters will be used. The difference be- 
tween the observed and expected redshift distributions can 
be attributed to the influence of omitted factors. Similar 
analysis can be repeated with richer data base and/or with 
available numerical simulations. 

2.3.1 Absorption in wall-like elements 

The first term in (2.2), n m , describes absorption associated 
with wall-like elements. Now these elements accumulate ~ 
50% of galaxies and their possible contribution cannot be 
rejected a priori. The observed mean separation of walls 
is about 40 - 60 ft _1 Mpc and the comoving linear number 
density of such elements is a w (0) «2.5 - 1.5T0 _2 /iMpc _1 
(LCRS1). These elements are formed during late evolution- 
ary stages and can be observed, primarily, at small redshifts, 
z <1 - 1.5. 

Taking into account expressions (2.3) & (2.9) we will 
approximate the mean dimensionless linear number density 
of absorbers associated with such elements by the two pa- 
rameter function 

n w {z) « k w E(c w ,z), (2.13) 

ccr w (p)a w _,, . , cd-2, x ,m 

k w = , E(c w , z) = exp[-c„(B (z) - 1)J, 

where the covering factor a w < 1, and the parameter c w 
characterize probability of formation of an absorption line 
in a wall and the period of wall formation, respectively. The 
observed wall-like structure elements are strongly disrupted 
(see, e.g., Fig 5 in Ramella et al. 1992) and therefore we can 
expect that a w -Cl, n w (0) = k w w 60a™. 

2.3.2 Absorption in filamentary elements 

The second term in (2.2), n/, describes absorption in fila- 
ments formed by DM, galaxies and intergalactic gas. The 
filamentary component of the structure accumulates « 50% 
of DM and galaxies and forms a joint network of structure 
between wall-like elements. As was noted above, the mean 
observed 2D surface density of galaxy filaments is estimated 
as <T/(0) « 0.01/i 2 Mp C - 2 . 

The mean free-path between filaments depends on their 
surface density, cry, and their diameters. To find it we con- 
sider the intersection of cylindrical filaments with typical 



radius of gaseous halo Rf and a random cylindrical core 
with radius r and length L. The mean number of such in- 
tersections can be found with standard methods (Kendall & 
Moran 1963, Buryak et al. 1994) as follows: 

(N int ) = -Ko-fL(r + R f ). 

The mean free-path between such cylindrical filaments along 
a line of sight can be approximated by 

(//) » L/<JV in( )|r=o = (na f R f )-\ (2.14) 

while the mean length of a line of sight within a filament is 
« 2R f . 

Taking into account expressions (2.10) we will approx- 
imate the mean dimensionless linear number density of ab- 
sorbers associated with filaments by a two parameter func- 
tion 

n s (z) » n f (l + z)B- 2 (z)E{c f , z) (2.15) 
«/ = -n-^Reffo-f(0), E(c f ,z) =cxp[-c / (B- 2 ( 2 )-l)], 

R °" = a ^ = Sff - w ^ (ooiife^) ^ lkpc - 

Here {Rf) is the mean radius of gaseous halo of a filament, 
the covering factor a/ (z) < 1 characterizes the probability of 
formation of an absorption line within a separate filament, 
and the factors 1 + z, and B~ 2 (z), describe the expected 
variation of surface density of absorbers associated with fila- 
ments due to the general expansion of the universe, merging, 
and expansion or compression of filaments. The parameter 
Cf characterizes the period of formation of the main fraction 
of observed filaments. 

The mean radius of gaseous halo of filaments, {Rf), can 
vary with the redshift and along the filament, because the 
observed galactic density varies along the filament as well. 
This means that our statistical estimate of averaged {Rf) 
is weighted by the matter distribution along filaments at 
different redshifts. 

More detailed models of absorbers associated with DM 
filaments can be developed. They depend however on many 
parameters which cannot be adequately determined from the 
available database. Hence, in this paper we will restrict our 
consideration to the two parametric model discussed above. 

2.3.3 Absorption in low mass structure elements 

The third term in (2.2), n pan , describes absorption by low 
mass structure elements. The available information about 
properties of this population at small redshifts is very lim- 
ited and we can estimate only the mean linear number den- 
sity of such pancakes n pan (0) = K, pa n. Taking into account 
expressions (2.12), we will approximate the mean dimen- 
sionless linear density of such absorbers, as 

n pa „(z) « K pan (l + z) 2 B~ 3 (z) E (c pan , z) , (2.16) 

E(c P an,z) = exp[-c pan (_B~ 2 (;z) - 1)]. 

Here the factors (1 + z) 2 , and B~ 3 {z), describe the expected 
variation of the density n pan (z), caused by the general ex- 
pansion, merging of pancakes, and their compression, and/or 
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expansion in transversal directions. The parameter K pan is 
the effective dimensionless free-path between such absorbers 
at z = 0, and c pan characterizes the period of formation of 
the main fraction of observed pancakes with N H i > N^j. 

The expressions (2.13), (2.15) & (2.16) give probable 
fitting relations to the observed evolution of absorbers asso- 
ciated with various types of structure elements. They allow 
us to estimate, in principle, the effective radius of gaseous 
component of filaments and the parameters of cosmological 
model Q. m and J1a- All fitting parameters have clear inter- 
pretation. 

Here we do not consider the possible contribution of any 
high density clouds, because they are usually embedded in 
DM filaments and pancakes. This contribution is small, at 
least at redshifts z < 3, and cannot be reliably singled out 
with available database, but it probably is more important 
at higher redshifts. 

2.4 Observed linear number density of structure 
elements 

The small scale clustering of absorption lines has been re- 
ported in many papers (see, e.g., Cristiani et al. 1995, 1996; 
Ulmer 1996), and it is especially important for metal sys- 
tems. It implies that a complex of nearby lines can be gen- 
erated in the same structure element. As we are interested 
in statistics of structure elements this factor must be taken 
into account. To do this we will use a technique developed 
earlier for the core-sampling method (Buryak et al. 1994; 
LCRS1), which allows us to select a Poison-like subsamplc 
of points from a general sample. The redshift dependence 
of the mean number density of the subsample characterizes 
the redshift evolution of uncorrelated absorbers, which can 
be identified with separate structure elements. 

This method uses separation of neighboring lines rather 
than the lines redshift itself, what attenuates somewhat the 
influence of selection effects inherent in individual spectra. 
If the number of lines in a sample is > 30 - 50 we can 
estimate the parameters of Poisson distribution with a rea- 
sonable precision 10 - 20%. Hence, we can use the avail- 
able catalogues, what is important for the metal systems, as 
number of such systems in individual absorption spectra is 
usually small. 

The major points of the method can be summarized as 
follows: 

(i) Separation between lines is characterized by the di- 
mensionless comoving distance 

Al = H Az/H(z). (2.17) 

(ii) A subsample of lines, in an interval z—dz < z < z+dz, 
taken from all spectra under investigation is organized into 
an 'equivalent single field' by combining the line separations 
Al one after the other along a line. 

(iii) Distribution of absorbers obtained in that way is as- 
sumed to be Poissonian for larger separations. ID cluster 
analysis is used to discriminate the Poisson-like subsample 
among the sample of points, and to find the number, and the 
mean linear number density of such points. The theoretical 
groundwork for such an approach was developed by Buryak, 
Dcmiariski, & Doroshkevich (1991). 

Using the standard ID clustering analysis, the mean 



number of Poisson points in the sample, Np, and their mean 
linear number density, np, are found by the maximum like- 
lihood fit to the relation 

\n(N l )=\n{Np)-n P R l , (2.18) 

where Ri and iV; are the variable linking length and the 
related number of clusters, respectively. Note that for a truly 
Poisson sample parameters NpEznp are related to the length 
of the 'equivalent single field', Do, defined by the first and 
the farthest point by 

Do/Np = n P . (2.19) 

Thus, the difference between the values Np and n p obtained 
from equation (2.18) and from equation (2.19), and varia- 
tions of the actual number of clusters along a straight line 
(2.18) give us a measure of error made due to the difference 
of the actual from the assumed (Poisson) distribution for 
absorbers along the line of sight. To decrease this error we 
use an automatic procedure, which finds the optimal range 
of Ri for the fit to equation (2.18). Both upper and lower 
limits of Ri were varied to obtain the best fit for the linear 
number density of absorbers n a bs{z). 

As a rule, the precision of this method increases with the 
number of lines, Nu ne (z), in the interval z — dz < z < z + dz. 
For Nn ne (z) > 50, any random variation does not exceed 
10%, and the real error is defined by the systematic varia- 
tions of the sample under investigation. For smaller Ni ine (z), 
random variations become essential. This condition restricts 
the choice of optimal interval to dz « 0.2 - 0.25. 



3 THE DATABASE 

The present analysis is based on spectra available in the 
literature. The list of such objects is given in Table 1. 

The distribution of lines over redshift is nonhomoge- 
neous and majority of lines are concentrated at z 3. Us- 
ing the technique described above, we can obtain the linear 
density of Poisson subsample of absorbers, at z < 3.2, with 
a reasonable precision of about of 10 - 20%. Distribution of 
absorbers, at z > 3.2, is based primarily on the 430 lines 
of QSO 0000-260 (Lu et al. 1996), and here the statistics 
of lines is rather poor. Inclusion of the spectrum of QSO 
1033-033 extends the redshift interval up to z ~ 4.4, but 
it cannot improve inadequate representativity of the sample 
at z > 3.2. Samples of poorer lines, with log Nhi < 13, are 
incomplete at all redshifts. 

A few samples of observed metal systems were analyzed. 
Here we present results obtained for the richest recent sam- 
ple (Vanden Berk et al. 1999). This catalogue contains 901 
separations between lines of neighboring metal systems from 
237 spectra of quasars at redshifts < z < 3.5. 

The published Ly-a lines with different Nhi were or- 
ganized into six samples listed in Table 2. All samples were 
supplemented by the HST data (Bahcall et al. 1993, 1996; 
Jannuzi et al. 1998). The samples with Nhi > 13.8, are 
associated mainly with the filamentary component of the 
structure, whereas samples with Nhi > 13, and Nhi > 12, 
are associated predominantly with the poor pancakes. 
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Table 1. QSO spectra from the literature 



Table 3. Fit parameters for the filamentary component at z <3. 
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1. Kulkarni et al. (1996), 2. Carswell et al. (1991), 3. Khare et al. 
(1997), 4. Kirkman & Tytlcr (1997), 5. Hu et al., (1995), 6. Lu et 
al. (1996), 7. Giallongo et al. (1993), 8. Cristiani et al. (1995), 9. 
Rodriguez et al. (1995), 10. Rauch et al. (1993), 11. Williger et 
al. (1994), 



Table 2. Samples of absorbers 
sample 

Qlt 
Q\f 
Q\l 
Q\l 

Q\l 

4 STATISTICAL ANALYSIS OF STRUCTURE 
EVOLUTION 

In this section the main results are presented for the redshift 
distribution of both Ly-a lines and metal systems. They con- 
firm that the redshift distribution of Ly-a lines is a super- 
position of several populations, with different evolutionary 
histories (Bahcall et al. 1996). Using the method and the 
fitting relations discussed in Sec. 2 for the analysis of sam- 
ples with different Nhi, we can roughly discriminate two 
populations of absorbers, and describe their evolution. Our 
estimates use essentially the HST data for z < 1.5. 

For three cosmological models the distribution of ab- 
sorbers was fitted to the expression 



Kf (1 + z) „ . , K 

n abs = n J E{c f , z) + - 



»(l + *) s 



E(c P 



B 2 (z) ~^ J '~' B3(z 

E(c,z)=e X p[~c-(B- 2 (z)-l)] 



(4.1) 



The function B(z) was introduced in (2.3). The main re- 
sults are plotted in Figs. 1-5, and the best fit parameters, 
Rjfi Cf , ftpun, and 

are listed in Tables 3-5. For samples of richer absorbers 
c P an=0, whereas for samples of poorer absorbers c pan ^ 
0.05 ±0.05 was found with large errors and, so, it is omitted 
in Table 4. 



i2 m 


n A 


nmet 
R eff 

h kpc 


c f 


h x kpc 


1.0 





40 ± 1.5 


0.07 ±0.03 


98 ±2.4 


0.5 





52 ± 1.6 


0.15 ±0.04 


118 ±3.3 


0.3 


0.7 


38 ± 1.7 


0.13 ±0.04 


69 ± 3.5 



4.1 The redshift distribution of metal system and 
stronger HI lines 

It can be expected that the redshift distribution of metal 
systems and stronger Ly-a lines is connected with the evo- 
lution of richer structure elements, which are seen in galaxy 
surveys as filaments and walls. To examine this hypothesis 
the function n abs (z) was found for the sample of metal sys- 
tems, and the samples of strongest Ly-a lines Q\\, and Q\l 8 . 
The analysis shows that in all cases no more than 10 - 15% 
of all lines can be attributed to the massive wall-like ele- 
ments, what is consistent with the large separation of walls 
(> 50/i -1 Mpc) and their strong disruption into a system of 
high density clouds, what decreases their covering factor. 

For all considered cosmological models the redshift dis- 
tributions of Ly-a absorbers, for samples Q\\, and Q}| 8 , are 
well fitted, at z < 3, to the one parameter function (4.1) with 
c f = 0, and = 0. The redshift distribution of metal sys- 
tems is well fitted to the two parameter function (4.1) with 
K pa „ = 0. The best parameters of these fits are listed in 
Table 3. They show that metal systems are predominantly 
formed within inner regions of gaseous halos of filaments, 
whose size, RTffi is about half the size of hydrogen halos, 
Rf/f- Both sizes are consistent with those observed directly 
at small z: Rhi ~ 100 - 150/i _1 kpc (Lanzetta et al. 1995), 
and R m et ~ 40 - 50ft _1 kpc (see, e.g., Le Brune et al. 1996), 
and to the expected size of DM halo (Bahcall et al. 1996). 
The precision of estimates of R e f / depends on the represen- 
tativity of the sample used. 

These results suggest that the redshift distribution, and 
the column density of absorbers are strongly correlated, and 
corroborate the identification of strong HI absorbers, and 
metal systems, with the filamentary component of the struc- 
ture of the universe. They show that such filaments were 
probably formed at redshifts z > 2, but their enrichment by 
metals could occur later, at redshifts z ft* 2 and smaller. 

At redshifts z > 3, the situation becomes more compli- 
cated, and the observed absorber distribution is no longer 
described by the one parameter fit plotted in Fig. 2 by dotted 
lines. To obtain a reasonable description, even for the red- 
shift distribution of stronger absorbers, at z > 3, the three 
parameter fit (4.1) with c pan = must be used. This indi- 
cates that at z > 3 progressively increasing part of stronger 
absorbers must be assigned to pancakes. 

The best parameters of these fits are listed in Table 
4, and fitting functions n a t s (z) are plotted in Fig. 2 by 
solid and long dashed lines. These estimates of c/ and R e ff 
crucially depend on the region of redshifts 2.5 < z < 3.3, 
and could be sensitive to the possible evolution of filaments. 
They show that the formation of such filaments described by 
the exponential term in (4.1) could occur at redshifts z > 1, 
while later the surface density of filaments progressively de- 
creases. The effective radii listed in Table 4 exceed by about 
2-3 times those listed in Table 3 for smaller redshifts. 
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Figure 1. n a (, s vs. z for metal systems for cosmological models 
with Q m = 0.3, Q A = 0.7 (a), Q m = 0.5, tt A = (b), and 

n m = i (c). 



This fact can be attributed to the progressive compression, 
small scale disruption, and relaxation of dark matter within 
and around filaments, which is accompanied by the progres- 
sive dissipation of gaseous halo of filaments. The observed 
evolution can also be partly assigned to the variations of 
background UV radiation with redshift. 

At redshifts z > 2 - 2.5, more and more essential frac- 
tion of strong absorbers is associated with rich pancakes, 
and, at z > 3.5, the pancakes become dominant, even for 
this population. This fact probably indicates the period of 
more rapid transformation of rich pancakes into filaments. 

This conclusion is based on poor statistics of absorbers, 
at z > 3, and should be tested with more representative 
sample of absorbers and with simulations. 

For our samples the parameters of standard fit to 



dN/dz = jV (l + z) 



(4.3) 



6.8 ± 0.5, 7z " 



12.7 ± 2, 7* 



0.35 ±0.25, 



1.5 ±0.3, 



(4.4) 



with a weak dependence on the considered cosmological 
models. 



Figure 2. The redshift distribution of absorbers for samples Q\\ 
(points & solid lines) and Q 14 s (triangles & long dashed lines), 
for cosmological models with f2 m = 0.3, Qa = 0.7 (a), Qm = 
0.5, SI a = (b), and fi m = 1 (d). Dot lines show the best two 
parameters fit for the filamentary component only at z <3. 



4.2 Redshift distribution of weaker Ly-a lines 

The samples of weaker absorbers are rich enough but, even 
so, the main fraction of observed absorbers is concentrated 
near z ~ 2.5 - 3.3. In Figs. 3 & 4 the functions n a t, s (z) 
are plotted, for the samples Q{%, Q\~2, and Q\\, to- 
gether with the four parameters fit (4.1). Both samples of 
poor absorbers with log Nhi < 13 are incomplete. The best 
fit parameters k/, c/, K pan , and c pan listed in Table 4 are 
weakly sensitive to the used cosmological models. 

These results show that, at redshifts z < 1.5 - 2, an es- 
sential fraction of weaker absorbers can be associated with 
the periphery of filaments and/or with possible membranes 
and bridges between branches of filaments, what is consis- 
tent with the expectations of Fernandes-Soto et al. (1996). 
But at high redshifts, the main fraction of poorer absorbers 
is certainly identified with the population of low mass pan- 
cakes and, for Nhi < 10 13 cm~ 2 , with absorbers formed 
within expanded regions. The small value of c van ~ 0.03 - 
0.05 suggests that these absorbers were formed at z m 4 - 5. 
These estimates are sensitive to the distribution of absorbers 
at z > 3, where the observed samples are not sufficiently 
representative. For these populations the random overlap- 
ping of poor absorbers in the redshift space, discussed by 
McGill (1990) and Levshakov & Kegel (1998), formation of 
absorbers within "minivoids" (Zhang et al. 1998; Dave et 
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Figure 3. The redshift distribution of absorbers for the sample 
Ql4 (points & solid lines) and the sample Q\ 2 (triangles & long 
dashed lines), for cosmological models with Q m = 0.3, SIa = 0.7 
(a), n m = 0.5, fl A = (b), and Q m = 1 (c). 

al. 1999) as well as the variations of UV background can 
distort the observed mass and temperature distribution of 
absorbers. As was noted in Sec. 2.1.2, at these redshifts, the 
more complicated three parameters function (2.11) provides 
better fit of observed distribution of absorbers. 

The essential variations of measured linear density, n a b s , 
with respect to the smooth fitting functions are seen in Figs. 
3 & 4 at z « 2.5 - 3. They can be caused, in part, by the 
poor statistics of absorbers, at z w 2.5, and z > 3.5. If the 
variations seen at z « 3 are real they can be connected 
with similar results obtained for the population of stronger 
absorbers and, thus, can indicate the period of fast transfor- 
mation of pancakes into filaments. They can also be caused 
by influence of local factors and, first of all, by variations of 
background UV radiation. 

The parameters of standard fit (4.3) are: 

N% an = 13.8 ± 0.6, 7 r" = 2.4 ± 0.2, for Q{ 3 2 , (4.5) 

N^ a " = 11.8 ± 0.6, 7 r n = 2.75 ± 0.2, for Q\%, (4.6) 

with a weak dependence on the parameters of considered 
cosmological models. The power index 7f an is close to the 
value found earlier (see, e.g., Carswell 1995; Cristiani et al. 
1996). 




z 



Figure 4. The redshift distribution of absorbers for the sample 
Ql4 (points & solid lines) and the sample Q]J (triangles & long 
dashed lines), for cosmological models with Q m = 0.3, Q A = 0.7 
(a), Q m = 0.5, Q A = (b), and Q m = 1 (c). 

4.3 Redshift distribution of Ly-a lines with b > 
40km/s and b < 20km/s 

The same technique can be applied to single out other sub- 
populations of absorbers such as subpopulations with lower 
Doppler parameter b, b < 20km/s, and larger, 6 > 40km/s. 
The main results of this analysis are plotted in Fig. 5, and 
the best fitting parameters are listed in Table 5. Unfor- 
tunately, they cannot be complemented by HST data at 
smaller redshifts. 

These subpopulations are not so representative and at 
z < 4 they contain only 640 lines with b < 20km/s and 727 
lines with b > 40km/s. This is the main source of the large 
errors of fitting parameters listed in Table 5. Nonetheless, 
these results demonstrate that, if the subpopulation of hot 
absorbers with b > 40km/s can be associated with pancakes, 
then the redshift distribution of cold absorbers with b < 
20km/s is similar to that found for filamentary component 
of the structure. 



5 SUMMARY AND DISCUSSION. 

In this paper the observed redshift distributions of Ly-a 
lines and metal systems are analyzed and interpreted in the 
framework of theoretical model of DM structure evolution 
(DD99). The observed evolution of absorbers depends on 
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Table 4. Fitting parameters of redshifts distribution of absorbers 
for three cosmological models. 



sample 



K f 



lOO/i-ikpc 



c f 





— 


1, n A =o 






Qli 


22 ± 1.3 


2.2 ±0.1 


0.34 ±0.14 


0.75 ± 0.12 


Qlf 


21 ± 1.2 


2.1 ±0.1 


0.33 ±0.15 


0.91 ± 0.11 


Qli 


43 ± 1.5 


4.3 ±0.1 


1.0 ±0.20 


4.2 ±0.11 




43 ± 1.3 


4.3 ±0.1 


1.0 ±0.20 


3.7±0.11 


Qli 


45 ± 1.6 


4.4 ±0.1 


1.1 ±0.22 


4.8 ±0.11 


o 12 


49 ±2.1 


4.8 ±0.2 


1.1 ±0.21 


4.0 ±0.09 





rim — 


0.5, Q A = 






Qli 


30 ± 1.5 


2.9 ±0.2 


0.66 ±0.21 


1.2 ±0.11 


Qlf 


23 ± 1.3 


2.3 ±0.1 


0.54 ±0.22 


1.6 ±0.14 


Qlt 


39 ± 1.5 


3.8 ±0.2 


1.82 ±0.37 


6.4 ±0.12 


Q\l 
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42 ±2.1 


4.2 ±0.2 


2.2 ±0.35 


4.4 ±0.11 



Table 5. Fit parameters for the subsamples of absorbers Q 
with b > 40km/s and b < 20km/s 





n A 


K f 




10^ ■ Cp an 






b > 40km/s 






1.0 





0.0 


2.4 ±2.4 


5.4 ±2.2 


0.5 





0.0 


1.9 ±2.0 


3.6 ±6.0 


0.3 


0.7 


0.0 


2.4 ±3.1 


7.7 ±4.1 






b < 20km/s 






1.0 





11 ± 4.2 


0.12 ±0.1 


1.0 ±0.5 


0.5 





16 ± 10 


0.0 


0.0 


0.3 


0.7 


10 ± 4.6 


0.0 


0.0 



several random factors and does not trace directly the evo- 
lution of the DM component. Our results show however that 
the observed redshift distributions of absorbers can be rea- 
sonably described even by the discussed model. 

The main results of our analysis can be summarized as 
follows: 

(i) The different redshift distributions of stronger and 
weaker absorbers favor the existence of two distinct subpop- 
ulations of Ly-a absorbers: a rapidly evolving subpopulation 
of weaker absorbers, that dominates at high redshifts, and 
more slowly evolving subpopulation of stronger absorbers, 
that dominates at low redshifts (Bahcall et al. 1996). This 
result coincides with theoretical expectations. 

(ii) The expected evolution of filamentary component of 
structure describes adequately the redshift distribution of 
metal systems, and stronger Ly-a lines, at z < 3. This strong 
correlation of evolutionary rate and the column density of 
absorbers suggests possible identification of stronger ab- 
sorbers with the filamentary component of observed galaxy 
distribution. 




Figure 5. The redshift distribution of 640 lines with b < 20km/s, 
(points, solid lines), and 727 lines with b > 40km/s for the sample 

:0.3,f7 A =0.7, (a), 



: 0.5, fl A = 0, (b), and Q m = 1, (c). 



(iii) The subpopulation of colder absorbers, with b < 
20km/s, can be probably related with the filamentary com- 
ponent of the structure. 

(iv) The rapid evolution of the linear number density of 
weaker absorbers, n a b s , can be naturally explained by the 
model of separate DM confined pancakes that are merging 
and expand or contract in transversal direction, as it is de- 
scribed by the theoretical model discussed above. At z < 2 
- 2.5, disruption of the compressed matter can essentially 
accelerate the evolution of this subpopulation. 

(v) The observed redshift distribution of absorbers de- 
pends on the influence of irregular local factors, such as 
variations of background UV radiation, and shock heating 
of gas produced by activity of galaxies. This impact is seen 
as irregular variations of the observed linear number density 
of absorbers around the smooth fitting curves. 

(vi) The obtained estimates of c/ can be interpreted as a 
probable pick of the rate of formation of filaments at z ~ 
2-2.5 whereas their enrichment by the metals can occur 
at z ~ 2 and smaller. The small values of c pan , found for 
all subpopulations under investigation, show that the main 
fraction of observed poor absorbers can be associated with 
DM structure elements formed at redshifts z > 5. 

The main quantitative characteristics of redshift distri- 
bution of absorbers, discussed above, admit natural interpre- 
tation and coincide with published estimates (see, e.g., Hu et 
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al. 1995; Lanzetta et al. 1995; Cristiani 1995, 1996; Le Brune 
1996; Bahcall 1996). They are roughly consistent with the 
expected evolution of DM structure and can be, in principle, 
used to test and to discriminate different cosmological mod- 
els. Now the best fit is obtained for the model with fl m ~ 
0.3, Qa=0.7, but the more representative database at red- 
shifts z ~ 2 - 2.5 and z > 3 is required for more reliable 
discrimination of cosmological models. 

On the other hand, the discussed rough model of struc- 
ture evolution ignores many details, such as, for example, 
the continual variation of morphological characteristics of 
structure elements, and does not describe the disruption of 
structure. This model could be improved, but then it is in- 
evitable to increase the number of fit parameters and, so, it 
will be justified only when richer sets of observed absorbers, 
especially at smaller and higher redshifts become available. 
This comment applies also to the more detailed description 
of expected variations of UV radiation and possible contri- 
bution of artificial caustics. Application of this approach to 
available simulations could test and improve it as well. 

The wall-like structure elements and large under dense 
regions found in distribution of galaxies at small redshifts 
cannot yet be reliably identified with the available database. 
The large scale modulation of redshift distribution of Ly-a 
lines found by Cristiani et al. (1997) and strong nonhomo- 
geneities, found recently at z < 2, by Williger et al. (1996), 
Quashnock et al. (1996, 1998), and Connolly et al. (1997) 
could be probably attributed to such rare extremely rich 
structure elements. Appearance of such elements is not for- 
bidden at any redshifts, but it can be expected that their 
number will rapidly decreases for larger redshifts. 
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